solve之所以慢,是因为它在默认情况下将浮子转换为理性主义,然后在这里试图求解一个非常大的多项式方程。如果您通过了rational=False,那么它会更快,但不会给出解决方案: print(solve(RR_equ, T, simplify=True, rational=False))
...
NotImplementedError: multiple generators [exp(2788.51/(T - 52.36)), exp(3346.65/(T - 57.84))]
No algorithms are implemented to solve equation -1.0 + 0 + 429133787.110874*exp(-2788.51/(T - 52.36)) + 782687755.282661*exp(-3346.65/(T - 57.84))solve函数用于寻找解析解,但您的方程不太可能有解析解: In [4]: RR_equ
Out[4]:
-2788.51 -3346.65
───────── ─────────
T - 52.36 T - 57.84
-1.0 + 429133787.110874⋅ℯ + 782687755.282661⋅ℯ = 0您想要使用的可能是类似于nsolve的东西,但是在这种情况下这是行不通的,因为方程由于非常大的指数项是病态的。通过快速查看lhs日志的图表,它实际上并不像这个特殊的方程有任何真正的解决方案(假设这就是你想要的)。
|